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Abstract 

A method is developed to derive simple relations among the reduced matrix 
elements of the quadrupole operator between low-lying collective states. As 
an example, the fourth order scalars of Q are considered. The accuracy and 
validity of the proposed relations is checked for the ECQF Hamiltonian of 
the IBM-1 in the whole parameter space of the Casten triangle. Furthermore 
these relations are successfully tested for low-lying collective states in nuclei 
for which all relevant data is available. 
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Microscopic shell model wave functions of collective nuclear states need a huge con- 
figurational space. However, experimental data indicates that there are comparably simple 
relations between the wave functions of different collective states, including the ground state. 
The wave functions of some excited states can be described by actions of one-body opera- 
tors on the ground state wave function with good accuracy. In even-even nuclei, where the 
ground state is a + state, the first 2 + state is given by the quadrupole operator Q acting 
on the ground state. The generalization of this concept has been named the Q"Pnonon 
approach Jl|-|7| . In this approach one describes the low- lying collective positive parity states 
of even-even nuclei in the basis of multiple Q-phonon excitations of the ground state, 1 ) , 

\L+,n) = N( L > n \Q...Q)^\0+) . (1) 

n 

In the framework of the IBM-1 ||,|§ it has been shown over the whole parameter space of the 
ECQF Hamiltonian [rO| , [TT| ] that each of the wave vectors of the yrast states can be described 
by only one multiple Q-phonon configuration with good accuracy [l]j4],[|], which has recently 
been confirmed by microscopic calculations in [[3]. 

The Q-phonon approximation implies the existence of selection rules for the matrix 
elements of the quadrupole operator. Thus, one finds that E2 transitions between Q-phonon 
configurations, that differ by more than one Q-phonon, are weak compared to transitions 
between those configurations that differ by only one Q-phonon. During the last years much 
data on 7-soft nuclei has been collected, especially in the A=130 mass region, which support 



these selection rules, e.g. [0-0 ■ 



The Q-phonon structure of the low-lying collective states allows one to obtain quadrupole 
shape invariants |14HT7| from rather few data. As an example we consider fourth order scalars 
obtained by coupling the four quadrupole operators in different ways. One obtains several 
different expressions for the fourth order quadrupole shape invariants in terms of only a few 
E2 matrix elements. These expressions can be used to derive approximate values of various 
observables, e.g., for the quadrupole moment of the first 2 + state or the lifetime of the first 
excited + state, from more easily accessible nuclear data. Such information is desirable for 
nuclei where complete experimental information about low-lying states is not - or not yet - 
available, as e.g. nuclei which are produced using rare isotope beams. 

There are three possibilities to couple the quadrupole operators to obtain fourth order 
scalars, 



P = (otKQ ■ Q)(Q ■ QM) , (2) 



^ ] = (oi\ [QQ] W [QQ] W v "'|oJ>, (3) 



i( 2 ) \nr>] ( 2 ) 



(0) 



gi 4) = (o^i[[Qg] w [Qg] w j^|o+), (4) 

The notation [. . .] W abbreviates the tensor coupling of two operators to angular momentum 
L. These three scalars are proportional to each other according to Dobaczewski, Rohozihski 



(0) 



and Srebrny in [18], if the Q-operators commute. Then one obtains the relations 
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In the IBM-1 the Q-operators do not commute. The effect of the noncommutativity of the 
components of the quadrupole operators scales however with 1/N and is therefore neglected 
in first order. Below, we will check the accuracy of Eq. (|5|) in the framework of IBM-1. 
In order to do this we decompose the scalars into sums over reduced matrix elements, 



?f= £Ki|g||2t><2+||g||o+) 

i,j,k 



x(0;||Q||2+)<2 fc +||Q||0+> , (6) 

gf = ^£<0+||g||2+)<2+||g||2+> 

o it j tk 

x(2+||Q||2+)(2+||g||0+) , (7) 



^ 4) = iE(o^iiQH2^)(2, + iigii4;) 



x(4+||g||2+)(2+||g||0+) . (8) 

Using Eqs. ©-(H), the quantities (@)-(|) have been calculated gridwise - using the code 
Phint |[TP| - for N=10 bosons over the whole IBM-1 symmetry space spanned by the 
ECQF-Hamiltonian [0,0 



H 



ECQF 



:i-C) n d -^Q-Q 



(9) 



The ECQF Hamiltonian interpolates between the symmetry limits of the IBM-1 using two 
structural parameters, £ and \- Here, is the <i-boson number operator and N is the total 
boson number. The parameter a has no structural meaning as it sets an absolute energy 
scale, and Q is the CQF quadrupole operator, both in the Hamiltonian and the E2 transition 
operator, 

l/e B T{E2) = Q = s + d + d + s + x[d + d) (2) , (10) 

depending on the structural parameter x, with —y/7/2 < x < 0; 6b is the effective boson 
charge. The result of this calculation is a near proportionality of the qf\ and in 
accordance with Eq. ([5]). 

In view of the selection rules of the g-phonon scheme, the sums reduce drastically. 

In the first approximation the set of E2 matrix elements necessary for the calculation of q± 
{n = 0, 2, 4) reduces to the following matrix elements, 

<2+||Q||4+> — > i = 1 , (11) 
<2+||g||2+)^ = l,2, (12) 
(2^||g||0 ? + ) — i = l,2,3. (13) 

The first three + states are taken into account, because the 0^ 3 -eigenstates of the ECQF 
Hamiltonian are mixtures of two- and three-g-phonon + configurations. Of course, if there 
are low-lying non- collective + states, the ECQF 0^ 3 eigenstates may refer to higher lying 
physical states. We have introduced the short notation Oqq by means of 
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(0fo\\Q\\J) 2 = (ti\\Q\\J)' + <<#\\Q\\J) 2 - (14) 

By using only the matrix elements (|TT| ) - (|T3|) in ®-(U) we see that i n eac h sum a factor 
{0f\ \Q\ \2f) 2 appears, which may be dropped, since we are interested in the ratios. Eqs. 
(!)"(!) now become: 

*1 0) = <2+||q||o+> 2 + <2+||q||o+ q > 2 , (is) 

4 2) = ^((2 1 -||Q||2+) 2 + (2+||Q||2+) 2 ) , (16) 
4 4) = ^(2+||g||4+) 2 , (17) 

where we use £4 to distinguish these quantities from the exact values. These quantities 
are also approximately proportional to each other, like the quantities q±. For arbitrary 
values of the boson number N the £4 values are related by factors defined by 

AO) _ 1 ,(2) _ 1 (4) 
c 02 c 04 

which depend on the boson number and the dynamical symmetry character. To obtain the 
values of the q^, we consider the U(5), SU(3) and 0(6) dynamical symmetry limits of the 
IBM-1 at first. The ECQF quadrupole operator ( |TU]) is used, and one obtains 

f 7N - 2 U(5) 
= J N(2N + 3) SU{3) , (19) 
[ N(N + A) 0(6) 

-±=(2iV + X 2 -2) U(5) 

= { i^( 4Ar + 3 ) 2 5f/ ( 3 ) , (20) 
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_(iV_ i)(iV + 5) 0(6) 

|(iV-l) U(5) 

tf> = { ^(JV-l)(2JV + 5) S77(3) . (21) 

|(iV-l)(iV + 5) 0(6) 



Comparing Eq. (0) and Eqs. ([I9| ) - (|2ll) one obtains proportionality factors for N — ► 00 

2 6 
^rTH ' C ^ = 35' (22) 

in agreement with the factors of Eq. (^). These values hold also for finite N in the 0(6) 
and the SU(3) dynamical symmetry limits when one neglects 1/iV 2 terms. Only in the U(5) 
limit a 1/N dependence is left, causing a small deviation from the limiting values. The 
values of the parameters with finite N differ slightly from those with iV — > 00. Using 
Eqs. (|IB])-(^1|) we obtain improved relations in the dynamical symmetry limits including 
the values of q^ for finite N. For nuclei far from symmetries one can calculate the exact 
q^ using Eq. ([18]) and interpolating in the IBM-1. We have done such calculation for the 
IBM-1 using the ECQF-Hamiltonian ([]). Fig. [1] shows the values of the parameter 
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1 - 



0/ 



(0i) = (02),(04) 



(23) 



for N=10 bosons. Using the limiting values for N — > oo results in a systematical error 
below 10%. We note that some deviations arise from our use of only the state for the 
Oqq configuration in this calculation. 

From Eqs. (|15|)-(|i"8|) one obtains two relations for the quadrupole moment of the 2f 
state: 



+ \2 



<2+||Q||2^ + (2+||Q||2 2 + > 



^3) ^x\\Q\m 

v c 02 c 04 / y 



(24) 



and 



32tt 
•35" 



' „N r oo " 
°02 c 04 

v c 02 c 04 . 



• £?(£2;4+^2+; 
- B(E2;2+ -> 2+ 



(25) 



(2+||Q||2+r + (2+||Q||2+>< 



JV 
02 



■y-S-«2f||g||0^ + (2f||g||0S g >' 

' C 02 



(26) 



Ql 



32tt 
.35~ 



-,iV 



~ 3S. [5^2; 2+^0+ 



c 02 



B(S2;0jQ-2+)]-S(S2;2+-2i 



(27) 



where Eqs. (|2"4|) , Q2"5D and (B^),(E^), respectively, differ only in notation. In a first approx- 
imation with c£/c§?=l, and if we define B{E2;2\ -> 2f) = 1/5(2^||Q||2^) 2 , we can write 
expression (Elf) in an intuitively interesting way: 



5(£2; 2+ -> 2+) + E(£2; 2+ -> 2+) = B(E2; 4+ -> 2+) . (28) 

A relation similar to (^) for N — > oo has been obtained in ||20|| , but was derived in a much 
less transparent way and was expressed using a rather difficult notation. Rewriting Eq. (j2~7D 
we get a relation for B{E2\ Oqq — > 2+), and a second relation by inserting (p7D in (]25|) . 
Extending our previous definitions []17|| of quadrupole shape invariants, we define now 



not only K±, but Kf\ and K^' , depending on the coupling. We want to obtain values, 
which characterize the nucleus and do not depend strongly on the coupling scheme. Thus, 
with q2=(Q\\Q ■ Q\0i), we introduce 



K4) 



K 



(0) 



K. 



(2) 



9f 



(2) 



K 



(4) 



3_5 ^ 
6 



(4) 



(29) 



(30) 



(31) 



These quantities are all equal if the quadrupole operators commute. We note that in the 
large N limit of the IBM-1 the theoretical values for the 1Q are 1 for the SU(3) and 
the 0(6), and 1.4 for the U(5) dynamical symmetry limit, distinguishing between /3-rigid 
and vibrational nuclei, respectively. Applying the above results to the i<Q leads to an 



approximation formula for that has already been obtained for N — » oo in ||T6"| 



K, 



A second approximation is 



4 0) 



(0) 



7 £(£2;4f -> 2 



10 B(E2; 2f -> OJ 



1 / = ^-appr. 



(32) 



7_ 

To 



+ E(E2; 2f ^ 2f ) 
5(£2;2f ^Of) 



(33) 



Due to K.f G [1,1-4] it emerges from Eq. (j33|) that, e.g., in the transition from 0(6) to 
SU(3), where K^=\, the value of Q 2 2+ /B(E2; 2f -> Of) rises from zero to 10/7, while the 

value of B{E2;2t -> 2f )/£(£2; 2f -> Of) drops from 10/7 to zero. Thus, these ratios 
characterize nicely the change of structure. 

In order to compare the relations with experimental data we considered nuclei near 
dynamical symmetries, for which all needed data is available. This data comes mostly from 
Coulomb excitation experiments by D. Cline and co-workers |1^J2T|J2^J26|J29[| . In Tables HLfT 



the results are given. The Os and Pt nuclei are considered to be 7- soft P2H24|1 or transitional 
between 7-soft and axially deformed nuclei, which is indicated by the large values of the 
quadrupole moments of the 2f states. As examples for vibrational nuclei Cd and Pd nuclei 
are shown, and Gd and Dy nuclei for the axially deformed case. We used values from 
the appropriate dynamical symmetry. 

In Table I the relations ( 18 ) are tested with satisfactory overall agreement. Additionally, 
the values of K2 Ppr ' are given in Table |. 

Table [Tl] shows the experimental values of Q 2 + and B[E2\ Oqq — ► 2f ) for the chosen 

nuclei, compared to the values obtained by the relations. The values of Q 2 +, obtained 

from the relations ( pop and (|2~T|), agree with the experimental values within the errors in 
most cases. A high accuracy of data is necessary for significant results, especially for the 
vibrator-like and the 7-soft nuclei, for which the quadrupole moments become very small. 

As an example for discrepancies, we consider 188 0s for which the B(E2; Of — > 2f ) value is 
very small, as expected for an 0(6) nucleus. The two-Q-phonon content of this state should 
therefore be very small. However, the values from Eqs. (p5|),(p7|) in Table [II] may refer to a 
higher lying + state with larger two-Q-phonon contribution, for which the lifetime is not 
known. Thus, with the missing E2 strength in 188 Os, the value of Q 2 + is underestimated by 
relation (0), while Eq. (p^| ) describes the quadrupole moment well. 

One finds significant deviations for other nuclei, too. For example, in 194 Pt the large 
B(E2\ Of — > 2f ) value indicates a two-Q-phonon structure for the Of state, in contradiction 
with the 0(6) prediction. Thus, this transition has been included in the calculation of 
the B(E2; Oqq — > 2f ) value. The value of K^ ppv ' = 0.8 in 192 0s is considerably smaller 
than its minimally allowed value: 1. This may be due to the small experimental value of 
B(E2- 4f — > 2f ). Also i^4 Ppr ' for 108 Pd is unexpectedly small, which does not support the 
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vibrational character of this nucleus. In the Cd isotopes considered the measured Q 2 + are 
smaller than expected from the relations. 

To summarize, we propose a simple method to derive sets of relations between the exper- 
imentally observable reduced matrix elements of the quadrupole operator. This approach 
is based on the use of the quadrupole shape invariants, the selection rules of the Q-phonon 
scheme and the fact that corrections from noncommutativity of the components of the 
quadrupole moment operator in the IBM-1 are small. As an example of the general scheme, 
fourth order Q-invariants of the ground state are given. One can apply the scheme also to 
higher order invariants, e.g. q 5 or q 6 , or to invariants built on excited states. The accuracy 
of the derived relations is checked for finite boson number iV over the whole parameter space 
of the ECQF-IBM-1 Hamiltonian and is shown to be rather good. A satisfactory agreement 
between data and theoretical relations has been obtained in many cases, but some exceptions 
clearly need further study. 

For fruitful discussions we thank A. Dewald, A. Gelberg, J. Jolie, T. Otsuka and Yu.V. 
Palchikov. One of the authors (P.B.) wants to thank the Institute for Nuclear Theory at 
the University of Washington for its hospitality during the final stages of this work, and 
R.J. thanks the Cologne University for support. We thank K. Jessen for the careful reading 
of this paper. This work was supported by the Deutsche Forschungsgemeinschaft under 
Contract No. Br 799/10-1. 
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TABLES 



TABLE I. Values of tf\ l/c Q 2 N tf ) and l/c 4 N t { ^ for various nuclei. The three values should 
agree according to Eq. [l8| 



data taken 
from 



,(o) 



3 2 b 2 l 



1/4 

[e 2 b 2 



t {2) 
t 4 



[e 2 b 2 



,(4) 
r 4 



m 

m 

\m 

pi 
2 



2.84(7) 
2.52(3) 
2.36(6) 
2.12(3) 
1.56(12) 
1.34(6) 



2.79(56) 
2.72(48) 
1.97(40) 
2.20(30) 
1.96(12) 
1.53(25) 



3.06(17) 
2.82(8) 
2.28(16) 
1.84(6) 
1.54(5) 
1.56(9) 



mm 



0.76(7) 
0.92(11) 
0.65(5) 
0.60(3) 



0.86(10) 
1.10(13) 
0.37(6) 
0.53(8) 



0.83(9) 
0.86(9) 
0.76(7) 
0.77(5) 



[1JH 
IIS 

in 



4.67(13) 
5.03(15) 
5.25(5) 
5.57(8) 



4.58(23) 
5.01(25) 
5.36(26) 
5.16(100) 



4.66(13) 
5.20(14) 
5.20(13) 
5.12(27) 
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TABLE II. Comparison of the quadrupole moments of the 2+ state and the reduced transition 
strengths of the 0+>q — ► 2+ transition for various nuclei. 



[e 2 b 2 ] 
Eg. (H) 



1.97^ 



[e 2 b 2 ] 



exp. 



[e 2 b 2 i 



+10 
-38 
+19 
-32 
+24 



B(E2;02 



3 2 b 2 l 



Egs. (| 



0.30+j 
< 0.10 


0.08±| 
0.24(11) 



.2u2l 



exp. 



0.040i 2 e 
0.0061+1 
0.014+1 
0.004+^ 
0.100+J 
0.02(1) 



[e 2 b 2 ] 

Eg. © 



1.80 
1.14 
0.56 


0.26(9) 



+7 

21 
+15 
30 
6 
-19 



1.72 
0.90 
0.84 



0.20+f 
0.24(18) 



+8 
-8 
+9 



+5" 

6 
+4 
8 



1.56 
1.20 
0.78+J 
< 0.01 
0.10(8) 



+TT 
11 
+11 

-11 



< 0.35 



0.20 



+15 
-20 





0.08+j 2 
0.50±? 7 
0.21(26) 



+TT 

11 
+11 
17 



22 
+9 
-20 



0.281 
0.20 
0.43(6) 
0.31(4) 



0.30 
0.38 
0.14(3) 
0.13(6) 



0.23 



0.24+1° 
0.35(5) 
0.18(3) 



0.20 
0.11 

0.27(8) 
0.26(6) 



0.14 



2 
+2 
2 



0.16] 
0.16(5) 
0.090(5) 



0.23 
0.35 


0.02(9) 



3.79(11) 
4.19(11) 
4.20(10) 
4.09(22) 



3.72(15) 
4.04(16) 
4.33(17) 
4.12(81) 



3.79(15) 
4.05(18) 
4.23(14) 
4.46(15) 



< 0.18 
0.17(20) 

< 0.09 




n.a. 
n.a. 
n.a. 
n.a. 



< 0.18 

< 0.28 
0.11(26) 

< 0.60 
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FIGURES 



N = 10 




FIG. 1. The deviations of the factors Cq 2 and Cq 4 from the limiting values, calculated gridwise 
over the whole ECQF symmetry space for N = 10 bosons. 



12 



